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Abstract 

By studying the structm'e of infrared divergences in a toy propagator in the 
replica approach to the Ising spin glass below Tc, we suggest a possible cancellation 
mechanism which could decrease the degree of singularity in the loop expansion. 
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1 Introduction 



The mean field approximation, when fluctuations are not taken into account, predicts a 
finite critical temperature for the Ising spin glass In the replica approach replica 
symmetry is spontaneously broken at Tc in a hierarchical and continuous way, yielding an 
ultrametric organization of equilibrium states (see [|, |], p for a review). 

Corrections to mean field theory up to third order in e = 6 — D in the paramagnetic 
phase have been carried out in [§, while problems still exist in the computation of such 
corrections below Tc. The difficulties lie in the complexity of the replica approach, which 
leads to very complicated bare propagators, with severe infrared singularities (see |[). 

A strong effort has been made in studying this difficult problem and many interesting 
results have been obtained. Among them, we recall the analysis of the modes with zero 
mass (see H, |l^), the study of the explicit breaking of replica symmetry (see [0, [T^), 
the study of some strong renormalization effects (see |n|) and, most recently, the 
derivation of a powerful method to compute the bare propagators (see |l^), which 
provides a clearer derivation of the formulae first presented in [|^ p. Unfortunately, the 
one loop contribution to the propagator has not been fully computed and even the lower 
critical dimension remains unknown. 

Ignoring the possible renormalization effects, a naive prediction for the lower critical 
dimension would be 3. Indeed, disregarding the singularity, which appears only in the 
zero overlap propagator and is not present in a small magnetic field, the leading infrared 
singularity of the bare propagators is p~^. However, it is not clear how much can really 
be inferred from these strong singularities in the framework of replicas. In particular, 
questions arise as to their origin and their consequences for measurable quantities. In fact, 
other approaches suggest that the lower critical dimension could be less that 3, or at least 
very near to 3. These follow from a computation of the free energy increase due to an 
interface among different phases and from numerical simulations |T^, |T^, . 



The aim of this paper is to try to take a step towards a deeper understanding of the 
structure of these divergences. We focus on how the singularities are induced in the 
bare propagators and we suggest a possible mechanism for their cancellation. 

The paper is organized as follows. After a brief summary of the replica approach, 
which can be skipped by the expert reader, we study a toy propagator defined using a field 
that explicitly breaks replica symmetry (already studied in [jl2| using a different method). 
For this propagator we are able to derive a differential equation which shows how the 
singularity due to continuous replica symmetry breaking cancels in the infrared limit 
and we illustrate how this cancellation results in a well-defined limit within the theory of 
distributions. Then we show how, in the framework of distributions, this propagator can 
be used to calculate the infrared behavior of a toy four-point function at tree level. We 
conclude our analysis by studying the structure of leading singularities in the full theory 
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and we discuss the possible generalization of the cancellation mechanism. 



2 The replica approach 

The Ising spin glass, or E-A model |1|], is defined by the following Hamiltonian, 

= -J2 Jij^i^j -hY^Si, (1) 

where Sj = +1, —1 are the spin variables, i = 1, A^, and () stands for nearest neighbors. 
The independent quenched parameters Jij are chosen from a Gaussian distribution with 
zero average and variance J'^ = 1/N. In what follows, the magnetic field h is taken different 
from zero, the limit h being performed after the limit oo. 

The study of the equilibrium properties of the quenched problem can be performed in 
the replica approach 0]. In this approach an effective theory is obtained by averaging over 
the disorder, indicated by an overbar in the following. Rather than do this for the free 
energy, the replica approach averages over the disorder the partition function of n copies 
(replicas) of the original model, n being analytically continued to zero at the end. The 
effective theory is then symmetric with respect to permutations of the n replicas (replica 
symmetry), with a. n x n matrix Q {Qab = Qta and Qaa = 0) for the order parameter. 

At mean field level, a second order transition occurs at the de Almeida-Thouless line, 
which terminates at = 1 for h = 0. On this line, roughly speaking, several states start 
to contribute to the Gibbs measure and ergodicity is lost. These results can be recovered, 
near T^, through the expansion of the effective free energy density F in powers of the 
(small) order parameter 

(3F[Q] + log2 + ^ = W[Q] = - hm -{^TtQ' + ^TrQ^ + 1 Y^iQat)'), (2) 

where t = — T and Tr stands for trace. In the framework of the Parisi ansatz , the 
saddle point of Q is looked for in a particular subspace using a hierarchical procedure. 
In this subspace, in which Q can be expressed in terms of a function q{x) defined in the 
interval [0, 1], the functional W[Q] results in 

W[q] = f\lq\x) - kxq^x) + 3q\x) f' q{y)dy) + ^q\x))dx. (3) 

Below Tc stationarity with respect q{x) yields the solution 

q{x) = x/2 < X < xi 

q{x) = q{xi) a;i < x < 1, (4) 

where Xi is defined by 2t — xi + xl/2 = 0. Replica symmetry, which in this framework 
requires q' = 0, is spontaneously broken below T^. 
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Replica symmetry breaking (RSB) is related to the probability of measuring a given 
value q for the overlap < Si >a< Si >b between two states, a and b, which differ by a finite 
amount in free energy. It has been shown that such a probability can be computed through 
the order parameter q{x), 

P{q) = = 2e(gi - g) + (1 - 2gi)5(g - gi). (5) 

In the metric dictated by the overlap these states are organized ultrametrically: given 
three states at least two overlaps are equal, the third being greater than the other two. 
This organization allows the use of a tree to express the overlap between different states. 
By putting the states at the end of the branches of a tree, the overlap between the states 
can be represented by the distance between the top root and the level of the point where 
the branches coincide. 

Before going beyond mean field, let us recall that the fluctuations of the order parameter 
Q around the RSB saddle point are usually divided into three families: longitudinal (L), 
anomalous (A) and replicon (R) (see ||16[ for the most recent and exhaustive analysis). 

The longitudinal modes are by definition invariant under the action of the symmetry 
group which leaves invariant the ansatz of and therefore correspond to fluctuations in 
q{x). On the other hand, the anomalous and replicon modes break even this n replica 
permutation group and are parametrized in term of functions of two and three variables, 
as explained in detail in |]16|. The conditions imposed by the breaking of this residual 
symmetry turn out to be strong enough to determine completely the R eigenvalues but not 
the L-A ones, where one has to explicitly solve the integral eigenvalue equations. 

Zero modes of the fluctuations around the mean field saddle point appear in each family, 
the rest of the spectrum being positive. A remarkable result is that in the replicon sector, 
where one has a closed expression for the eigenvalues A(x, fc, /), one finds that zero modes 
are present for a finite range of eigenvectors, i.e. 

A(a;, x, x) = 0. (6) 

This can also be seen by explicit differentiation of the saddle point equation and using the 
fact that replica symmetry breaking is continuous (as first shown in |^). 

For future reference, let us also recall that an explicit RSB can be introduced in the 
theory (a deep analysis on the nature of the explicit RSB can be found in |rT|) by adding 
to the effective free energy the term 







1 

q{x)e{x)dx. (7) 



A finite conjugate field e induces a shift in the order parameter q{x) which provides a 
kind of infrared regulator because it induces a gap proportional to the slope of e in the 



spectrum (as shown in ||12|| ). 

In the next section we consider such a field e as an external source in order to perform 
a detailed analysis of the infrared limit in the replica approach. 
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3 The "projected" theory 

For the time being, let us define a propagator in the subspace identified by q{x). To define 
such propagator, we add to the functional W[q] a kinetic term and consider how a small 
external conjugate field e, explicitly breaking replica symmetry, perturbs the mean field 
solution. We have the following theory for a field 6q{x;p), in which x is a continuous 
internal degree of freedom, 

W[q + 6q] / {6q{x;p)Ydx + Sq{x;p)e{x;p)dx, (8) 

4 JO Jo 

and we introduce the bare propagator G{x,y;p) through the relation 

6q{x;p)=f G{x,y]p)e{y]p)dy. (9) 
Jo 

For small e the equation for q{x) leads the following equation (see 0), 

p^G{x, y;p) + 2 T qiz)G{z, y; p)dz + 2q{x) /' G{z, y; p)dz = 25{x - y). (10) 

Jo Jx 

After repeated differentiation of the replica variable, we obtain 

^, d^G{xj-p) _ 2^,^^)^^^^ y. ^ 2g"(x) £ G{z, y- p)dz = 26" {x - y). (11) 

In what follows we are mainly interested in the diagonal sector x = y < Xi when p ^ 0, 
so we focus on this case. For p = the previous equation leads to the solution 

G{x,y;0) = -26"{x-y), (12) 

while the solution to these equations for finite p, and x,y < Xi is 



26{x — y) e 



G{x,y;p) = — + g{x,y;p), (13) 

where the function g{x, y; p) is not singular in the limit p — (|^). 

We observe that propagator already computed has been already studied in |T^, using 
a different method, for the purpose of discussing the regularization induced by e. 

The results (0, [T^) are interesting for two reasons. On the one hand, (|I^) shows how 
the small momentum behavior of this propagator can be cast into the form of a distribution. 
The integral kernel of the zero momentum equation for G is only min{g(a;), q{y)} because 
the strictly diagonal contribution of the kernel vanishes on the mean field saddle point. 
This implies that eq. (0) is its inverse, which shows that a small e{x) > 0, independent 
of space, induces a shift of q{x) for x < Xi only through —2e{x)" > 0. The propagator 
induced by e for x < xi is massive. 
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On the other hand, (|13|) allows us to understand how continuous replica symmetry 
breaking gives rise before xi to the diagonal singularity for small p. In fact, when 
we add the kinetic term of order p"^ on the diagonal, the absence of a strictly diagonal 
contribution in the kernel implies a contribution on the diagonal of G. But now, to 
keep the off-diagonal elements of the product of the two matrices zero, a diverging off- 
diagonal contribution in G is also needed. Such a contribution can only be p~^ with an 
exponential prefactor because it has to be smooth and coalesces for p = 0. 

We learn that care is needed in the limit p 0. This limit should be considered 
in the sense of distributions, paying attention to the cross-over between x and p. A finite 
momentum induces a regularization of the distribution which appears through two singular 
terms, which cancel in the limit of small p. 

To proceed in this analysis let us now consider the propagator with a source e 7^ 0. 
We know from the previous analysis (at e = 0) that for small p there is an off-diagonal 
contribution of width p and order p~^ that cancels the p~^ singularity and leads to a massive 
propagator. The analysis for p = and small e is similar with e'/q' playing the role of p"^. 
At zero momentum the propagator does not diverge as e'~^ but is instead given by 

G(.,y-,O).= '-0^^-'^^. (14) 

6e{y) q'{x). 

Let us use this result to compute the four point function at tree level. By taking two 
derivatives with respect to the conjugate field e we obtain 

GH. „(, , 0, ^ ^^^1.^0 ^ _, n.-.)n.-.Or(.-.) . ,,,, 

de{z)de{w) q'{X)^ 

This result, derived in detail in 0, shows that in the four point function the infrared 
contributions from the two diagrams with four external legs, the one from the quartic 
vertex and the one from two cubic vertices with a propagator flowing between, are similar 
but do not cancel. 

The absence of the complete cancellation is not too surprising for this propagator. In 
fact this propagator is massive and expresses the response of the order parameter to an 
explicit breaking of the replica symmetry. It is similar to the longitudinal propagator in 
an 0{N) model, which is regular in the infrared limit. However it is remarkable that in 
the long wavelength limit the propagator forces the two diagrams with four external legs, 
which are very different in structure, towards the same kind of contribution. This is exactly 
what happens in an 0{N) model, where the cancellation of a different four point function 
(transverse) is required by a Ward identity. 

The possibility of using distributions to investigate the infrared limit of the complete 
propagators and to compute the complete four point function is appealing. Let us analyze 
the structure of the leading singularities of the full theory. 
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4 A preliminary analysis of the full theory 



Our aim here is to investigate the possibihty of extending the results obtained for the toy 
propagator, a two index object, to the complete propagators. 

We have performed a preliminary study of these propagators using the results in 0, § . 
The propagators for the full theory are defined through the inverse of the mass matrix with 
a diagonal kinetic term 

G.p,^sip) =< SQap{p)6Q,s{-p) >= + ^^^) ^ • (16) 

Because the replicon eigenvalues are known in a closed form the Green functions are usually 
split (we follow the notation introduced in |^) into two contributions, 

""GZ^ip) and "^^GZM, (17) 

where a n/3 = x if qa/s = q{x), ada = 1, and x = anP,y = 'ynS,zi= max{a 07, a fl 5} 

and Z2 = max{/3 fl 7, /3 fl 6}. 

We are interested in the divergences of order because we know from that the 
singularity is confined to strictly zero overlap and disappears in a small magnetic field. 

From the full propagators given in P] one obtains in the long wavelegth limit in the case 

of < x,y < xi (see ^ 



z < x,y 
u < z 

Zi, Z2> X 
Zi, Z2> X 

where the last two formulae have already been given in [0] and u = min{x,y}, v = 
max{x,y}. Let us analyze these results. The singularities are confined to the diag- 
onal L-A contribution, where the ultrametric prefactor is u~^. In the R propagator the 
singularities are of order and the ultrametric prefactor is 

We are tempted to understand all these singularities as being generated by the same 
mechanism as those of the toy propagator of the previous section. That is, the zero modes 
in the R fluctuations {X{x,x,x) = 0) induce the singularity in the diagonal (upper 
indices) L-A propagator. 

However, it seems that there is a difference with respect to the toy propagator. In fact 
in the toy propagator the singularity is canceled because of the diagonal p~^ singularity. 



^^G^^(p) 



e p 
up^ 

u — v 

e p 

up^ 
1 

xp^ 
1 

9 9 ' 
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a b c da b c da c b d 



Figure 1: From the left to the right we show the ultrametric trees which correspond to the 
propagators GH.G^l and G^f^j. 

opposite in sign. In the full theory the ultrametric prefactor of the diagonal R contribution 
is not the same as the L-A one, the sign also being the same. 

To proceed in this analysis, let us observe that when one has four ultrametric indices, 
there are several different possibilities of arranging them on an ultrametric tree. In fig. |I| 
we show the trees corresponding to three of the possibilities. The three graphs correspond 
to the three possibilities of having two equal indices in the four index propagator. In the 
notation used for the propagators, where at least two indices must be equal, the graphs 
correspond, from left to right, to ^^Glf /"^G%f and ^^Gff^j +^^1"^^^. When the two equal 
indices are the upper ones (x = y), there can be both an R and an L contribution to 
the propagator (third graph in fig. [l|), apart for the boundaries ^1,2:2 = ^ where the R 
contribution is zero. 

Because ^"^G^^^^ ~ does not depend on ,21,^2 > 2;, we are interested in the volume 
of the corresponding L-A diagonal subspace, i.e. the sum 

E ''^(^M (19) 

abed 

where 9 is equal to 1 only in the diagonal L-A region. This sum can be performed using 
the general result of Mezard in which a sum over replica indices can be transformed into a 



sum over all the possibilities of arranging the indices on an ultrametric tree . 

A sum over several replica variables is equivalent to a sum of different contributions, 
each one corresponding to a possibility of arranging the variables on an ultrametric tree. 
Each tree has a weight that is given by the number of possible permutation of the tree 
multiplied by a factor which depends on the structure of the tree. If we specialize this result 
to the case of four indices, we obtain that for each node with three branches the factor 
is X while for each node with four branches the factor is 2x^. The previous prescription 
holds for all different indices. When the indices coincide pairwise there is an additional 
multiplicative factor (—1)™, where m is the number of pairs. 

Using these prescriptions to compute (|19|) we obtain 

2x^ + 2(1 - xf + 4x(l - x) - 4a; - 4(1 - a;) + 2 = 0, (20) 



8 



i.e., the volume of the L-A diagonal subspace vanishes. The volume of the R subspace 
zi, Z2 > X is given by 

2(l-x)2-4(l-x) + 2 = 2a;^ (21) 
while the volume of the off-diagonal subspace is 

a; + 4(1 - x) - 4 + 2x = -x, (22) 

which is exactly the factor necessary for the L-A propagator to have the same x dependence 
found in the R one! 

Let us conclude by considering the projection of the four index propagators in a two 
index subspace. In other words we want to define a two index object, through the four- 
index replica propagators, by integrating the lower indices. This operation is defined by the 
sum of all the propagators Gab,cd over all replicas a, b, c, d at fixed x = aHb and y = cHd. 
Using distributions for small p we find that 

EW.) - -2x^fc^ + 2xf^i:^-^:^^i^) 

^ P \ xp^ X J 

- -26"{x-y) 

= G{x,y;0). (23) 



5 Conclusions 

In this paper we have analyzed a toy propagator that expresses how a small external field, 
explicitly breaking replica symmetry, induces a perturbation on the order parameter q{x). 
This propagator, defined in the subspace identified by g(x), turns out to be, as expected, 
the projection of the complete set of propagators in this subspace. 

We have shown that for finite p this propagator is essentially given by two contributions, 
their singularities canceling for p ~ 0. This cancellation results in a well-defined infrared 
limit within the theory of distributions that we have used to extract the infrared behavior 
of a toy four point function at tree level. 

We have considered some aspects of the full theory. In particular we have analyzed 
the ultrametric structure of the subspace where the fluctuations are of order p'^, which 
corresponds to the third graph in fig. 0. This subspace has a global volume equal to zero 
and includes, apart from the boundaries, the replicon subspace. Because the volume of the 
off-diagonal subspace is —x one might conjecture that, by casting the infrared behavior of 
the propagators within the theory of distributions, these singularities cancel. 

A careful analysis is necessary, and work is in progress in this direction. 

It is a pleasure for us to thank C. de Dominicis, I. Kondor and T. Temesvari for 
interesting discussions and R. Monasson and M. Virasoro for a useful collaboration. We 
are grateful also to D. Lancaster and M. Potters for the careful correction of the manuscript 
and for very useful suggestions. 
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A Projected propagator 

The equation of state at first order in e is 

which leads to following equation for G{x, y) 

p^Gix, y;p) + 2 T qiz)Giz, y; p)dz + 2q{x) f' Giz, y; p)dz = 25(x - y). (25) 



The solution of eq. ( ^Sf ) for finite p can be achieved by solving the inhomogeneous equations, 
adding the most general solution of the homogeneous equations and searching for the 
coefficients which solve eq. ( ^Sf ) . We obtain for x,y < Xi 



26{x - y) e-\--y\/P , 

G{x,y;p} = + g{x,y;p} (26) 

pi 

where the function g{x,y) is 

g{x, y; p) = ^++e+("+^)/P + ^~e-("+?')/P + ^+-(e-'"-^l/P + e+l'-^l/^) 



p- 


(1- 


a;i) 










1 p + 


p3 
(1- 


Xi) 


{P + 


(1 - Xi))e+(^-i/p) + (p- 

g+(3:^l/p) 


(l-,Ti))e- 


-(a-l/p) 


p- 


p3 
[1- 


Xi) 


{P + 


(1 - a;i))e+(^i/p) + (p - 

g-(a;i/p) 


(1 - xi))e- 


-(a;i/p) 



-- 

^ p3 (p+(i_a;^))e+(-i/p) + (p_(l_Xi))e-(-i/f)' ^ ^ 

as obtained in by using the longitudinal eigenvalues. As can be seen from eq. (^), 
when X or y goes beyond the breakpoint Xi the solution, apart the delta function, does not 
depend any more on this variable. For x and y > Xi we then obtain 

^ 25ix - y) _ 2_ e(-^/P) - e^-^/^) 

^^'^'^^ p2 p2(p+(i_3,^))e{xi/p) + (p_(i_2,^))e-(xi/p)- ^^^^ 



In the limit p — this gives 



, 25ix-y) 2 1 

G X, y- p) ^ ^-Y^ - --. r. 29 

p^ p^ (1 — Xi) 



10 



B Projected four point function 

The aim of this appendix is to show how the result for four point function can be derived 
through the "projected" propagator. By a derivative with respect to e we obtain 

dG{x, y; 0) ^ _ 1 d dq{x) 

de{z) ^\'{xydxde{z) 



and by an additional derivative we obtain 

G'SL(^,y,z,«;) = 5"{x-y){35"'{x-z)^ 



d dq{x) 



q'(x)^ dx de{w) 



q'[xY dx de{w) q'[xY dx"^ de(w 

S"{x - y) 



q'{xf 
5"{x - y) 



q'{xY 



{■i5"'{x - z)5"'{x -w) + 5"{x - z)5""{x - w)) 
{26"\x - z)6"\x - w)) (31) 



In the last step the contributions of the two diagrams with four external legs, which do not 
cancel because a multiplicity factor 3, can be recognized. 

This result can also be obtained using diagrams. The expansion of W around the mean 
field saddle point in the subspace identified by q[x) gives a cubic and a quartic vertex 

(3) 5{x- y)d{z -y) + 2 permut. 
V^>{x,y,z) = 

V^'Hx,y,z,t) = (32) 

Using these vertices, together with the projected propagator, to calculate the four point 
amputated function at the tree level and zero momenta we find 

G%{x, y, z, t) = -4!y(^) {x, y, z, t) + 3(3!)^ J^' dudvV^'^ {x, y, u)G{u, v; 0)V^''> {v, z, t) 
= 5{x - y)S{y - z)S{z - t) (-2 + 6) 7^ (33) 

which is exactly the same result obtained by differentiation on the propagator after cutting 
of the external legs. 
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C Infrared limit of the full propagators 



In this section we use the results given in for the full Gaussian propagators as the 
starting point to explicitly derive their behaviour in the infrared limit. Using the results 
and the notation presented in P] the propagators are given as integrals of a kernel F^'" 
weighted with an ultrametric measure. 

In the L-A sector we obtain that for small p the formulae for the leading contribution to 
the kernels in |^ simplify. We consider first the sector k < u < v. The leading contribution 
of order p~^ is independent of k and disappears in the integral by the derivative. However 
when k is within order p of u the kernel becomes of order In fact we have 

-Ti''^r^--'^i-- + -J (34) 
kdk " k p^ ^ p k' ^ ' 

which, for the propagator, leads to 

2z — u — v 

'^G7. ^ (35) 



For u < k < V we find that the kernel for small p does not contribute to leading behavior 
of the second and third propagator, which are then given only by the kernel k < u < v 



-\^-y\ 



e p 
up 

The last L-A propagator (the diagonal) is then 

1 

up 

In the R sector the divergent contribution to the propagator for < x < a;i is given by 

Rr.xx ^ r ^ dk2_d^ 4 

21^2 / u / i^h A^2 I U2 I U2 



"^^^Z. ^ (37) 



ki Jx k2 dkidk2 4p2 + kj + k^ - 2x^ 

2x{zi — x) r2x{z2—x) 



pAxyz-i — x) rzx(Z2—X) 

/ drii / dr]2- 

Jo Jo I 



Jo Jo (4j92 + V1+V2) 



j^2p2 



{31 



that is independent of zi, Z2, if they are greater than a; by a finite amount. 
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